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Abstract We present multi-wavelength imagery along with the X-ray emission maps of 
merger remnant galaxy NGC 1316 to study the content of dust and its association with other 
phases of ISM. Color-index maps as well as extinction maps derived for this galaxy reveals an 
intricate and complex morphology of dust i.e., in the inner part it exist in the form of a promi- 
nent lane while at about 6-7 kpc it apparently forms an arc-like pattern extended along the 
North-East direction. Apart from this, several other clumps and knots are also evident in these 
maps. Dust emission mapped using Spitzer detection at 8 fjm indicate even complex structure 
its morphology. Extinction curve derived over optical through near-IR bands exhibit identi- 
cal extinction properties of dust and can be assessed from the parallely running extinction 
curve with that for the Galaxy. Quantum of dust estimated from optical extinction measure- 
ment is found to be 2.13 x lO^M© , while that from IRAS flux densities is 2.11 x lO'^M© 
and from integrated flux densities at 24 pm, 70/im and 160 pm from MIPS is 3.2 X IO^Mq, 
significantly larger than the estimates from the optical extinction. High resolution Chandra 
observations of this merger remnant system have provided with an unprecedented look at the 
complex nature of distribution of X-ray emission that closely matches with that of the ion- 
ized gas and to some extent with the dust distribution also. X-ray color-color plot of the 80 
resolved X-ray sources within optical D25 extent of the galaxy enabled us to separate them in 
different classes. 

Key words: galaxies: individual (NGC 1316), galaxies: ISM, X-rays: ISM, (ISM:) dust, 
extinction, galaxies: elliptical and lenticular, cD 
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1 INTRODUCTION 

Multi-wavelength data on early-type galaxies from both ground based as well as space borne telescopes 
have helped us to shed some light on the origin of multiphase ISM in this class of objects (Goudfrooij et al., 
1994; Rampazzo et al., 2005; Patil et al., 2007; Cappellari et al., 2011; Finkelman et al., 2010b, 2012). The 
general picture that emerges from the past studies is that either the cold/warm ISM observed in these galax- 
ies could be of external origin, accreted through merging and/or close-encounter episodes with neighboring 
galaxies, or be the by-product of cooling of the hot ISM that has originated from internal stellar mass loss, 
e.g., cooling flows (Kim & Fabbiano, 2003). Systematic study of these galaxies selected from different re- 
gions, like fields, groups and clusters, etc. delineate that the relative contribution of the two competitive 
process i.e., external versus internal, do not follow a general rule regarding their origin but varies among 
the objects. Though much attention has been paid to probe the properties and origin of ISM in early-type 
galaxies with the spectroscopic study of the dynamics and chemical abundance of the gas and stars (Bertola, 
1987; Kormendy & Djorgovski, 1989; Sarzi et al., 2006; Annibali et al., 2010), however, due to the limited 
spectral coverage of these investigations and limited knowledge of the interplay between different phases 
of ISM, it was not possible to arrive at a decisive conclusion regarding the true nature of the ISM. 

Detailed investigation of morphological forms and orientation of dust and ISM in these galaxies could 
help us in getting insight of the process of their origin. The unusual optical signatures discernible in the 
form of dust-lanes, shells, tidal features, double nuclei, etc. are the direct evidences indicating towards 
the merger like event that the host might have experienced in the past (Bertola, 1987; Goudfrooij et al., 
2004; Patil et al., 2007). If it is having internal origin through the evolution of single stellar population, 
instead of getting oriented in the form of disks, lanes, etc., it will be evenly distributed throughout the 
galaxy (Finkelman et al., 2012). Another indication comes from the fact that the amount of dust content of 
such galaxies quantified using the optical extinction measurements as well as from IRAS flux densities is 
always larger by several factors than that expected from the mass loss of evolved stars (Goudfrooij et al., 
2004; Dewangan et al., 1999; Patil et al., 2007; Finkelman et al., 2010b, 2012). Furthermore, kinematical 
studies of this class of galaxies have revealed that the motion and orientation of the gas in many systems 
is decoupled from the stellar rotation, and hence provides an additional evidence for their external origin 
(Bertola, 1987; Kormendy & Djorgovski, 1989; Caon et al., 2001). Thus, multi-frequency study of various 
components of ISM in E/SO galaxies with dust lanes are important to investigate the nature and origin of 
the ISM in this class of galaxies. 

NGC 1316 is a peculiar SO galaxy, with numerous tidal tails and shells and many pronounced dust 
patches including a dust lane oriented along its optical minor axis. In addition to the this, NGC 1316 
also hosts Ha filaments, strong shells (Malin et al., 1983), many filaments and loops (Schweizer & Seitzer, 
1988). All these features strongly suggest that NGC 1316 have experienced a strong merger event 
(Terlevich & Forbes, 2002). In X-rays NGC 1316 has been observed with Einstein (Fabbiano et al., 1992), 
ROSAT PSPC (Feigelson et al., 1995) and with CHANDRA (Kim & Fabbiano, 2003). X-ray study shows 
presence of a low luminosity AGN with Lx ^ 5.0 x 10'^^ erg s^^ in energy range 0.3-8.0 keV 
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Fig. 1 B-band image of NGC 1316, overlaid on which are the Chandra X-ray point sources 
detected within ACIS-S3 chip. This figure clearly reveals a set of shells, ripples as well as tails; 
signatures of merger remnant 

(Kim & Fabbiano, 2003). In addition to the nuclear activity, this galaxy also exhibit extended X ray emis- 
sion with many substructures. NGC 1316 is one among the nearest giant radio galaxy with well defined 
core-jet-lobe structure. The lobes show a very smooth distribution of light with sharp boundaries. The 
bridge of emission is found to be significantly displaced from the centre of NGC 1316, again an evidence 
for the strong merger event (Ekers et al., 1983). Goudfrooij et al. (2001) based on the study of the bright 
globular clusters have estimated age of this merger remnant to be ^3 Gyr Thus, NGC 1 3 1 6 is a good target 
to investigate dust extinction properties and to compare those with the canonical grains in the MiUcy Way. 

In this paper we discuss the results of multicolour CCD imaging that involve broad band B V R I 
and narrow band Hq, study. We mainly focus on the dust distribution, extinction properties and compare 
morphology of ionized gas with that of dust in NGC 1316. We also discuss the results of reanalysis of 
archival Chandra observations of NGC 1316 that involve imaging as well as spectroscopy of target galaxy. 
This paper is organized as follows. In section 2 we describe the optical and X-ray observations and the 
data analysis. Section 3 discuss the dust extinction properties and morphology of ionized gas. The results 
of spatial and spectral analysis along with those of discrete sources are also discussed in this section. In 
section 4 we discuss the issue of origin of dust as well as its association with other phases of ISM. Finally, 
we summarize our results in section 5. We adapt a distance of D=25Mpc throughout this analysis. 
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Table 1 Global parameters of NGC 1316. 



Parameter 



Value 



Alternate names 

RADEC 
Morphology 
MagBT 
D25 

Distance(Mpc) 
Redshift(z) 
Radial velocity(km/s) 
Effective radius(kpc) 
IRAS flux density(Jy) 

1 .4GHz flux density(mjy) 



Fornax A.IRAS 03208-3723,PGC 12651 

ES0357-22,NVSS J032241-371225 

03:22:41.7; -37:12:30 

SABO(s)pec 

9.42 

12'.0x8'5 
23.6 

0.005871 
1653±12 
7.08 

0.330±042(12 ^im);0.43±0.02(24 fim) 
3.07±0.03(60^im);8.11±1.99(100 fim) 
254.7 




Fig.2 Central 2'. 5 x 2'. 5 region of NGC 1316, (left panel) shows (B-V) colour index map, 
darker shades in which delineates the dust extinguished regions. Cross indicates the position of 
the galaxy center recognized by NED, while {right panel) image shows the B-band extinction 
map of NGC 1316. 
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2 OBSERVATIONS AND DATA PREPARATION 

2.1 Optical data and near-IR data 

Deep, high S/N CCD image on NGC 1316 in B, V, R, I, Ha filters presented in this paper were ac- 
quired from the NED and were originally observed with the Cerro Tololo Inter-American Observatory 
(C770)/University of Michigan Curtis Schmidt 0.6/0.9 m telescope (see Mackie & Fabbiano (1998) for de- 
tails). The detector was a Thomson 1024 x 1024 CCD, each pixel of 19/im with a spatial resolution of 
l."83. These are processed images i.e., bias-subtracted, flat-fielded, co-added after registering centroids of 
bright field stars and exposure corrected (Mackie & Fabbiano, 1998). Sky background was estimated using 
the box method and was then subtracted from the respective passband image. These geometrically aligned, 
background subtracted images were then used for the dust extinction study in the target galaxy. One of such 
image in B passband is shown in Figure 1 . This figure reveals several shells, ridges and ripples around the 
main galaxy at larger radii, indicative of the merger signatures. With an objective to extend the study of dust 
extinction in near-IR bands, we have also acquired J, H & Ks band images on NGC 1316 from archive of 
Two Micron AU Sky Survey (2MASS) observatory. 

2.2 X ray data 

Though X-ray emission properties of NGC 1316 have already been reported using the observations from 
Chandra (Kim & Fabbiano, 2003; Lanz et al., 2010), XMM-Newton (Isobe et al., 2006; Lanz et al., 2010), 
SUZAKU (Konami et al., 2010), Einstein (Fabbiano et al., 1992), as wefl as ROSAT (Kim et al., 1998), how- 
ever, with an objective to examine association of hot gas with dust and ionized gas in this merger remnant 
SO galaxy NGC 1316, we have made use of high resolution X-ray data on NGC 1316 from the archive 
of Chandra observatory. NGC 1316 was observed by Chandra (Obs. ID 2022) on 2001-04-17 with the 
ACIS-S3 as an aimpoint for an effective exposure of 30.0 ks. 

Standard tasks available within Chandra Interactive Analysis of Observations (CIAO version 4.2.0) 
and recent calibration files provided by the Chandra X-ray Centre (CXC) (CALDB version 4.3.0) were 
employed for the analysis of X-ray data. These data sets were first filtered for the periods of high background 
emission using 3cr chpping of the 0.3 - lO.OkeV light curve extracted from the chip with a binning of 260 s. 
This resulted in to the net exposure time of 25.0 ks. For background subtraction we used properly scaled 
blank sky background files provided by the CXC. Point sources recorded on chip S3 were detected using the 
wavdetect tool within CIAO adopting a detection threshold of 10^^ and scale parameter covering 6 steps 
between 1 to 32 pixels. This has enabled us to detect a total of 86 discrete sources within chip S3 out of 
which 80 sources were found to lie in the optical D25 region of NGC 1316 and their distribution is shown 
in figure 1 . 
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3 RESULTS 

3.1 Dust properties 

3.1.1 Dust Extinction 

Though NGC 1316 is known to host dust features since their first time detection by (Schweizer, 1980), 
quantitative analysis of extinction properties of dust in this galaxy are not available in the literature. The 
amount of extinction and its effect on observed colors strongly depend on the size distribution, structure 
and chemical composition of the dust grains. For knowing the extent and spatial distribution of dust in 
target galaxy, we generated color index maps (B-V), (B-R), (B-1), etc. for NGC 1316 by comparing light 
distribution in the geometrically aligned, seeing matched broad band B,V,R,I images. Figure 2 {left panel) 
shows the (B-V) color-index map of NGC 1316, where darker shades represent the dust occupied regions 
and are consistent with those found by Schweizer (1980). This figure reveals an internal prominent dust 
lane along the optical minor axis of NGC 1316, which then forms an arc like pattern at about 6-7 kpc. In 
addition to these main features, several other filaments and dust clumps are also evident in this figure. 

Comparison of light distribution in extinguished part of early-type galaxies with that in the absence 
of dust allows us to investigate extinction properties of dust and its wavelength dependent nature (Brosch, 
1987). This is possible by deriving dust free smooth models of target galaxies in different passbands. Here, 
dust free models of target galaxy NGC 1316 were generated by fitting ellipses to the isophotes in optical 
broadband images using the ISOPHOTE package within IRAF (see Patil et al. 2007; Vagshette et al. 2012a 
for details). The position angle, ellipticity and center coordinates were kept free during this fit till the signal 
reaches 3(7 of the background. Before fitting the ellipses, regions occupying dust features and foreground 
stars, evident in colour index maps, were masked and ignored during the fit. Dust free models thus generated 
were then used to quantify the wavelength dependent nature of dust extinction, extinction curve, using the 
relation 



where Aa is the amount of total extinction in a particular pass band (B,V,R,I) measured in magnitude scale, 
while Ia,o6s and l\,modei are the observed (attenuated) and undistinguished light intensities in a given pass 
band, respectively. One of such extinction maps derived for NGC 1316 is shown in Figure 2 (right panel). 
This figure confirms the unusual, intriguing and clumpy morphology of dust in the target galaxy, as was 
evident in the (B-V) colour index map. In the inner part, dust appears like a well-defined lane oriented along 
optical minor axis and then at about 6-7 kpc it appear in the form of an arc in the North-East direction. In 
addition to this, several knots, dust patches and clumps are also evident in this figure. 

To investigate quantitative properties of dust extinction and to examine its wavelength dependent na- 
ture, total extinction values were measured over the dust-extinguished region in each of the band by sliding 
a 5 x5 box on the dust occupied regions. The numerical values of local extinctions in different pass bands 
(Aa) were then used to derive the extinction values in each pass band by fitting a linear regression between 
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Fig. 3 Extinction curve for NGC 13 16 

the total extinction A\ and the selective extinction E{B — V) ^ Ab — Ay - The best-fitting slopes of these 
regressions along with their associated uncertainties were subsequently used to derive the and hence the 
extinction curve. The extinction curve derived over optical-to-near-lR region of the electromagnetic spec- 
trum for the dust occupied regions of the Fornax A (NGC 1316) galaxy is shown in Figure 3. From this figure 
it is clear that the dust extinction curve plotted between the measured values of Ra ^= e(b-v) ) versus in- 
verse of wavelength of their measurement varies linearly, as was observed for the Galactic extinction curve 
(Mathis et al., 1977) and in extinction curves derived for different galaxies derived by (Goudfrooij et al., 
1994; Patil et al., 2001, 2003, 2007, 2009; Finkelman et al., 2008, 2010b,a, 2012; Vagshette et al., 2012a). 
This extinction curve for NGC 1316 imply that, though their extinction properties are identical with the 
Galactic as well as extragalactic environment, the grain particles responsible for the extinction of optical 
and near-IR light are relatively smaller than the canonical grains in the Milky Way with their average grain 
size equal to < ^J-^ >= 0.93 ± 0.08 with average Ry value equal to 2.91±0.12. 

3.1.2 Dust mass estimation 

Total dust extinction measured above were used to quantify the dust content of the NGC 1316. For this 
we integrated the dust column density S^; over the area (A) occupied by the dust in the extinction map. 
Assuming that chemical composition of extragalactic dust grains is uniform throughout the galaxy and is 
similar to that of the dust in our Galaxy, one can compute total dust content of a galaxy using available mod- 
els for the chemical composition and shape distribution of the dust grains. Several grain models, invoking 
different size distribution functions and mixtures of chemical compositions, are proposed in this context. As 
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our dust extinction study is restricted to the optical and near-IR bands only, therefore it is possible to em- 
ploy a simple model that can account for the Galactic extinction curve. For this we adopt the (Mathis et al., 
1977) two-component model comprising of spherical silicate and graphite grains with an adequate mixture 
of grains, that successfully reproduce the Milky Way extinction in the optical-UV part of the electromag- 
netic spectrum. This model assumes uncoated refractory particles with power-law size distribution. While 
the details of the dust estimation are described in Patil et al. (2007), the dust mass can be estimated by 
integrating the dust column density over the dust occupied area (A) as 



graphite grains, the hydrogen number density; is the dust column density; is the specific grain 
size density (^3 g cm^'^); Id represents the dust column length along the line of sight; and amax and Omm 
represent the upper and lower cutoffs of the size distribution, respectively. 

Using the measured values of average extinction in V band (< Ay >), we estimate the dust content of 
NGC 1316 to be equal to 2.13 x 10^ Mq. As this method assumes the screening effect of dust, therefore 
is insensitive to the component of dust that is diffusely distributed throughout the galaxy, therefore, the 
clumpy dust masses calculated above gives only a rough estimate of the dust content of target galaxies. The 
uncertainties involved due to the lower and upper cutoffs of the grain size may also worsen this estimate. 

The dust content can alternatively be estimated adopting the far-IR emission flux densities measured at 
60/U,m and lOO/^tm using the IRAS satellite. We estimated the dust grain temperature using the IRAS flux 
densities at 60/im and lOO/^tm using the relation Td = 49 (y-§^^ K and is found to be equal to 32.5 K. 
The dust mass estimated from these IRAS flux densities is found to be equal to 2.11 x 10^ M0 and is an 
order of magnitude higher than that estimated using optical extinction method. This discrepancy in the two 
estimates is due to the fact that, optical method is insensitive to the intermix component of dust, while IRAS 
can record this component. We have also employed the integrated MIPS flux densities at 24/im, 70yum and 
160/xm to quantify the dust content of this galaxy and is found to be equal to 3.21 x 10^ Mq. This estimation 
is perfect agreement with that reported by Lanz et al. (2010) and Draine et al. (2007). 

3.2 Morphology of warm gas 

With an objective to examine the association of dust with ionized gas, we have derived H„ + [N II] emission 
map of NGC 1316 following the method discussed by Mackie & Fabbiano (1998). This was done by sub- 
tracting the properly scaled, sky-subtracted continuum image from the sky-subtracted emission-line image 
Mackie & Fabbiano (1998). The adopted scaling was calculated from a linear least-square fit to residual of 
the field stars. Spatial distribution of the ionized gas within NGC 1316 is shown in Figure 4(/e/r panel) that 
closely matches with that of the dust evident in (B-V) colour map as well as X-ray emission map in 0.3-3.0 
keV band smoothed with 2 sigma gaussian (Figure 4). The arc of ionized gas in this figure appears more 
prominently. Figure 4{right panel) indicates the spatial correspondence of dust with the ionized gas. 




where n{a)da nijAia 



3.5 



da), with a being the grain size, Ai the overall abundance of silicate and 
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Fig. 4 (left panel): Ha emission map with 3a smoothing, contours overlaid on which (white) 
represent the dust morphology, (right panel): 0.3-3.0 keV X-ray image smoothed with 3(7 over- 
laid with Ha emission contours (pink) and dust contours (white) 

3.3 X-ray emission 

X-ray emission map is the best tool to study the morphology hot gas distribution in a galaxy. With an 
objective to examine the association of dust and Ha emitting ionized gas with the X-ray emitting hot 
gas, we have derived 0.3 -3.0 keV X-ray emission map from the Chandra observations and is shown in 
Figure 4 {right panel). This figure represents the background subtracted, exposure corrected, point source 
removed, adaptively 3cr smoothed 0.3-3.0keV Chandra image of NGC 1316. X-ray emission map derived 
for NGC 1316 is in perfect agreement with that derived by Kim et al. (1998) using ROSAT data and that 
of (Kim & Fabbiano, 2003; Lanz et al., 2010) using Chandra data. Figure 4 clearly reveals morphological 
similarity between the distribution of ionized gas and the hot gas exhibiting the strong correspondence be- 
tween the two. For comparison we overlay the contours of Ha emission map on the X-ray image. Both these 
morphologies closely follow that of the dust (Figure 4 right panel) in this merger remnant galaxy pointing 
towards common origin of all the three phases of ISM i.e., through the merger event. X-ray emission in the 
energy band 0.3-3.0 keV appears in more extended form compared to that of the ionized gas and shows a 
very disturbed structure along with filaments and patchy halos, perhaps X-ray cavities. This figure exhibit 
that orientation of the dust lane in inner region is slightly off of the hot and warm gas. 

As was discussed above, we have detected 80 discrete X-ray sources within the optical D25 region of this 
galaxy. Presence of these sources are always thought to be linked with the star formation history of their 
host in general and formation scenario of the host galaxy in particular (Irwin et al., 2003; Colbert et al., 
2004; Vagshette et al., 2012b). X-ray color plot for such sources has been taken as an efficient tool for 
investigating characteristics of individual source (Prestwich et al., 2003; Fabbiano, 2006). Position of in- 
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= (H-M)/(Sh-Mh-H). 



dividual source in this X-ray color plot clearly delineate their intrinsic nature and hence an unique tool 
to classify them in different classes, perhaps linked to their formation scenario Vagshette et al. (2012b). 
Therefore, with an objective to investigate their X-ray characteristics, we have derived X-ray color plot for 
these sources. For this, we extracted background subtracted X-ray photons from individual source in three 
different bands namely, soft (S, 0.3-1.0 keV), medium (M, 1.0-2.0keV) and hard (H, 2.0-lO.OkeV), using 
the task dmextmct available within CIAO. Then hardness ratio of individual sources was estimated using 
H21 = s+mIh ^31 = s+m^-Ih Vagshette et al. (2012b). The plot between the X-ray hard (H31) 
versus soft color (H21) is shown in Figure 5. From this figure it is apparent that, majority of the sources 
resolved in this merger remnant galaxy are of normal low mass X-ray binary (LMXB) type and are shown 
by green circle (Vagshette et al., 2012b). Four of them are of high mass X-ray binary type (HMXB), four of 
heavily absorbed type, one of super-soft type and one is of heavily absorbed very hard source, perhaps the 
nuclear AGN candidate. In the case of star forming galaxies, generally, population of super soft sources due 
to the accretion of matter on to a white dwarf compact companion is observed to be more. These detection 
and classifications are in perfect agreement with those reported by Kim & Fabbiano (2003). Their approach 
for this classification is different. 

With an objective to examine global properties of hot gas distribution in the target galaxy, we have 
extracted a point source removed, background subtracted X-ray spectrum of the photons from within op- 
tical D25 region of NGC 1316. Central 20"region was excluded from this to avoid the contribution from 
AGN. This spectrum was fitted following the standard statistics within XSPEC version 12.6.0q, with 
an absorbed single temperature model (MEKAL). However, the fit exhibited residuals particularly in the 
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Fig. 6 Azimuthally averaged projected radial profile of (a) Temperature, (b) abundance, (c) 
electron density, (d) pressure measured with Chandra. 

higher energy range. Therefore, we added a power law component to it, to constrain the emission from 
the unresolved discrete sources. This resulted in to a relatively better fit with value close to 1.62 with 
83 dof. During this fit, the hydrogen column density was frozen at the Galactic value (2.40 x 10^° cm~^) 
(Dickey & Lockman, 1990). The best-fit resulted in to kTcoo; = 0.18±0.03, kT/,ot = 0.62±0.02 and photon 
index r=0.45±0.10. 

Given that the hot gas is spatially associated with the ionized gas as well as dust, with a view to examine 
the temperature structure of hot gas distributed within NGC 1316, we have performed spatially resolved 
spectral analysis of the X-ray photons. As was observed earlier, soft X-ray photons have extended distribu- 
tion relative to harder photons, therefore we extracted 0.3-3.0keV X-ray photons from 8 different annuli 
centered on the X-ray peak of NGC 1316. Width of each of the annulus was set such as to get roughly same 
number of counts to validate the statistics. The background spectrum was extracted from the exposure 
corrected blank sky frame provided by the CXC. Source spectra, background spectra, photon-weighted re- 
sponse files and photon-weighted effective area files were generated for each of the annulus using the CIAO 
task specextmct. Spectra extracted from each of the annulus were fitted with a single temperature thermal 
plasma model apec with neutral hydrogen column density fixed at the Galactic value (Dickey & Lockman, 
1990). Temperature, metal abundance and normalization were kept free during the spectral fitting. The 
resultant, azimuthally averaged projected radial temperature profile derived for NGC 1316 is shown in 
Figure 6. From this profile it is evident that the X-ray photons distributed within NGC 1316 shows tem- 
perature structure in the sense that its temperature is dependent on the radial distance of their distribution. 
It monotonically increases as a function of distance, similar to that observed in the case of cooling flow 
galaxies (Gastaldeflo et al., 2007; Rasmussen & Ponman, 2007; David et al., 2009; Pandge et al., 2012). A 
jump in the temperature profile is evident at about 1.'2 and is perhaps linked with the dust absorption. 
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4 DISCUSSION 

The issue of origin of dust in early-type galaxies is highly controversial. Its internal origin in early- 
type galaxies assumes main contribution from the atmospheres of asymptotic giant branch (AGB) stars. 
Supernovae (SNe) have also been recognized as a potential candidate for the injection of dust and gas in 
to the ISM, where they can supply ^ 0.1 to 1 MQof dust (Dwek et al., 2007). The dust injected in to the 
ISM at the same time is processed by the forward and reverse shocks in the hot gas swept up by SNe and 
hence may get destroyed (Nozawa et al., 2010). With the observed SNe rate in NGC 1316 and considering 
the two competitive processes of formation and simultaneous destruction of dust grains, we can estimate 
total content of dust accumulated by the galaxy over its life time by solving the equation, 

dMd{t) dMd^, 



dt dt 



where ^^g''^'-*-* is the dust accumulation rate; ^^^f^ is the rate of dust injection by the SNe and stellar winds; 
Md{t) is the available dust at a given time t; and r^^ is the rate of dust grain destruction. We estimate the 
Hfe time of dust grains of radius a against collision due to electrons, hot protons and a-particles adopting 
relation given by Draine & Salpeter (1979) and is found to be 4.1 x 10^^ yr^^. By assuming gas-to-dust 
ratio of 100 we estimate the dust injection rate of 0.09 M0 yr^^ and hence total build up (Figure 7) of dust 
mass in NGC 1316 over its life time of ^^3 Gyr to be equal to 2.0 x 10^ Mq. 

Comparison of this estimate with the true dust content of this galaxy derived by using IR luminosity, 
indicate that the injected amount of dust is shorter by an order of magnitude than the actual content. This dis- 
crepancy may enhance if we include estimate from the observations at the sub mm wavelength (Chini et al., 
1995). Our estimate based on the integrated MIPS fluxes estimate leads to a value of 3.21 x 10^ Mq and is 
significantly larger than the quantum acquired by the galaxy internally. Thus, the internal supply of dust by 
the mass loss from evolved stars is not sufficient to account for the observed amount of dust in this galaxy 
and hence supports the external origin of ISM through merger like events. Lanz et al. (2010) demonstrate 
that this galaxy has formed out of the merging of a late type galaxy with a supply of 2—4 x 10^ Mq gas, 
sufficiently enough to account for the observed values of dust as well as gas in this galaxy. There are several 
other evidences that indicate towards the formation of Fornax A through a merger like event. The first and 
most important evidence is the presence of shells around the main optical part of the galaxy. Lanz et al. 
(2010) further demonstrate that, like Centaurus A, Fornax A is also a remnant of a recent merger exhibiting 
a stronger nuclear activity due to the presence of a 10** MQblack hole at the centre. 

Systematic study of NGC 1316 has enabled us to examine the spatial correspondence between the 
interstellar dust. Ha emitting ionized gas as well as the X-ray emitting hot gas. In Figure 8 we illustrate the 
spatial association of all the phases of ISM, even by including dust emission at S/^m from the Spitzer data. 
Dust occupied regions mapped through the emission from PAH observed at 8/xm using Spitzer is found 
to coincide well with the X-ray emission in the central 7.5 kpc region. Further, the X-ray emission from 
this galaxy is co-spatial with the ionized gas and is having extension along its minor axis and is perhaps 
due to the hot plasma ejected out from the supernovae in the central region. Similar features are also seen 
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Fig. 7 Dust accumulation by NGC 1316 over its life time through the internal process of its 
formation 

in some other cases of star forming galaxies (Strickland et al., 2004; Vagshette et al., 2012a). Multiphase 
association of ISM in this merger remnant galaxy NGC 1316 is more clearly illustrated in figure 8. 

5 CONCLUSIONS 

We present broad-band as well as narrow-band images along with the X-ray emission maps of NGC 1316 
with an objective to study the content of dust and its association with warm and hot gas. Results derived 
from this study are: 

- Color-index maps as well as optical extinction maps derived for NGC 1316 reveals significantly different 
morphology of dust. It hosts a prominent dust lane oriented along its optical minor axis in the inner 
region, then at about 6 kpc it appears in the form of an arc along the North-East direction. Apart from 
this, several other dust knots/clumpy features are also evident in this galaxy. The dust emission maps 
derived from the Spitzer observations at 8 /im indicate even complex features. 

- Extinction curve derived for the range between optical {B, V, R!kl) to near-IR (J, HSzKs) bands exhibit 
identical physical properties of dust grains in the merger remnant NGC 1 3 1 6 as those of the Milky Way, 
and can be confirmed from the parallel running extinction curves. However, relatively smaller value of 
Rv indicate towards smaller grains responsible for the observed nature of dust extinction. 

- Quantum of dust measured from the optical extinction is about 2.13 x 10^ M0 and is significantly 
smaller than the estimate from IRAS flux densities as well as from integrated flux densities at 24/xm, 
70/im and 160/im from MIPS. This in turn point towards the external origin of the dust. 

- High resolution Chandra observations of this merger remnant Fornax A system have provided with an 
unprecedented look at the complex features of distribution of hot gas in this galaxy. 
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Fig.8 Association of multiphase ISM in NGC 1316; (a) Hydrocarbon emission at 8/im from 
Spitzer, (b) symmetric stellar light distribution 2MASS K-band image, (c) optical light distribu- 
tion in B-band, in which obscuration due to dust is also evidence, (d) NUV GALEX image, (e)is 
the continuum subtracted Ha emission map and (f) tri-color map with hydrocarbon emission at 
8/.im (red), optical light distribution (green) and diffuse X-ray emission (0.3-3.0 keV) (blue). All 
the images are in 4 x 4 arcmin^ region. 



- Morphologies of the hot gas and the ionized gas mapped through the X-ray emission and Ha emission 
maps, respectively, exhibit surprising similarity in the two and to some extent with the dust features 
also, hinting towards common origin of all the three components of ISM. 

- Spectrum of the diffuse X-ray emission within optical D25 region of the point sources removed image 
is well constrained by a two-component (apec + power-law) model, where the harder component due 
to the unresolved population of sources was represented by the power-law. 

- We could resolve 80 point sources, in contrast to 81 detected by Kim & Fabbiano (2003), within D25, 
majority of which are like normal LMXBs with 1.4 Mq neutron star accreting sources and are associ- 
ated with the galaxy. X-ray color-color plots for the resolved sources in target galaxy exhibit a significant 
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difference in nature of the X-ray emitting sources. We found one heavily absorbed harder source and 
only one supersoft source, consistent with the findings of Kim & Fabbiano (2003). 
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